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ABSTRACT 

We apply the lDelove fc Bildstenl l)2003|) isentropic models for donors in ultracompact low-mass X-ray 
binaries to the AM CVn population of ultracompact, interacting binaries. The mass-radius relations 
of these systems' donors in the mass range of interest (M 2 < Q.IMq) are not single- valued, but 
parameterized by the donor's specific entropy. This produces a range in the relationships between 
system observables, such as orbital period, P or b, and mass transfer rate, M. For a reasonable range 
in donor specific entropy, M can range over several orders of magnitude at fixed P or b- We determine 
the unique relation between M and M 2 in the AM CVn systems with known donor to accretor mass 
ratios, q = M 2 jM\. We use structural arguments, as well as each system's photometric behavior, to 
place limits on M and M 2 in each. Most systems allow a factor of about 3 variation in M, although 
V803 Cen, if the current estimates of its q are accurate, is an exception and must have M 2 ~ 0.02Mq 
and M « 1O _1O M0 yr . Our donor models also constrain each donor's core temperature, T c , range 
and correlate T n with M 2 . We exam ine how variations in donor specific entropy across the white 
dwarf family l)Nelemans et al.l|2001b|) of AM CVn systems affects this population's current galactic 
distribution. Allowing for donors that are not fully degenerate produces a shift in systems towards 
longer P rb and higher M increasing the parameter space in which these systems can be found. This 
shift increases the fraction of systems whose P or b is long enough that their gravity wave (GW) signal is 
obscured by the background of detached double white dwarf binaries that dominate the GW spectrum 
below a frequency w 2 mHz. 

Subject headings: binaries: close — gravitational waves — stars: AM CVn 



1. INTRODUCTION 

The AM CVn class of variable stars are He-rich ob- 
jects which show a striking absence of H and photo- 
metric and spe ctroscopic var iations with periods of « 
300-4000 s (see IWarnerl 119951 for an overview). These 
periods appear to be orbital and the commonly held 
model for these systems is that of a mass-transferring bi- 
nary with a low-mass He white dwarf (WD) donor and a 
C/O or He WD accretor l)Paczvriskll967t|Faulkner et alJ 
I1972J) . This model is supported by several lines of ev- 
idence: t he spectra are dominated by double peaked 
lines (e.g. , IMaxsh et abllTflflfll ICroot et al]l200l. imply- 
ing the presence of an accretion disk l|Nather et al.lfl98lL 
iNasser et al.ll200"l|) . rapid photometric flickerin g is ob- 
served, suggesting ongoing ma ss transfer (e.g. IWarnerl 
119721 IWarner fc Robinsonll972D . and, in several systems, 
the spectral line profiles vary with a stable, regular pe- 
riod coincident with one of the primary ph otometric pe- 
riods , indicating this is the s ystem's P or b l|Nather et al.l 
I198H iNelemans et aLll2001al) . These systems' large ob- 



serv ed spectral lin e bro adening requires a co mpact accre- 
tor llNather et alJll98ltlPatterson et al.lll993URuiz et alJ 
120011 ICroot et al.l 120011) while their weak X-ray emis- 
sion indicates the accretor is a WD l|Nather et al.lll98ll 
IRuiz et al.ll2lffllT:IIJllalll99l . 

Because AM CVn systems are undergoing mass trans- 
fer, we can constrain the donor's structure from the or- 
bital period, P or b- For stable mass transfer, the donor's 
radius, P2, must equal its Roche lobe's radius, Rl, which 
can be approximated as 
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v Mi + M 2 

when q = M 2 jM x < 0.8 l|Paczvhskilll967|) . Here M x and 
M2 are the masses of the accretor and donor respectively 
and a is the orbital separation (throughout, we assume 
that the orbit is circular) . Combining Rl = R 2 with Ke- 
pler's third law leads to the period-mean density relation, 
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In analyzing and modeling AM CVn systems, it is cus- 
tomary to specify a mass-radius, M-R, relation for the 
donor producing a one-to-one relation between M 2 and 

Porb- 
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Two such M -R relations have been utilized extensively, 
each motivated by a different AM CVn formation chan- 
nel. In both of these ch annels, several common envelope 
(CE) episodes (see, e.g.. lPaczvnskll976tlYungelson et alJ 
Il993|) are involved in producing a close binary pair con- 
sisting of a degenerate WD accrctor and a low-mass 
He core remnant. The latter may eventually becomes 
the donor when and if gravity wave angular momentum 
losses bring the pair into contact. In the first channel, 
which we will refer to as the WD-channel, this He ob- 
ject never ignited after the end of the CE phase, and the 
donor makes contact as a (possibly semi-) degenerate He 
WD. In this case, the donor is typi cally modeled using a 
T = M-R relation such as that oflZapolskv fc Salpete"rl 
l|1969fl (e.g.. INelemans et aUl2001b|) . In the other chan- 
nel, which we refer to as the He-star channel, the He 
object is able to ignite He and makes contact as a He- 
burning star dSavoniie et "aTlll986t iTutukov fc Fedoroval 
Il98?l lErgma fc Fedoroval 1 199ft INelemans et al.ll2001bfl . 
For this channel. |Nelemans et alJ (l2001bf) used a fit to the 
evolut ion track of model 1.1 from ITutukov fc Fedorova 
1989) to obtain a donor M-R relation, while | Warner 



1995 ) has also quoted a M-R fit to the iSavoniie et al 



, 1986) ev olutionary calculation. In additio n to these two 
channels. iPodsiadlowski et alJ <|2002tl20"o"3|) have recently 



proposed a scenario involving stable mass transfer from 
a H-depleted main sequence star that could also produce 
AM CVn systems. The potential contribution to the AM 
CVn population from this channel has not be systemat- 
ically explored nor have M-R relations for the donors 
produced in this channel been reported. 

These sets of M-R relations have been used exten- 
sively to a pproximate the properties of known sys- 
rWaxneHI1995HProvenca,l et a,UI1997t IMarsK 



terns [e.g.. IWarneiJ 11995} iProvencal et~aTlll997t iMarshl 
1999} iNasser et al.l 120011 lEl-Khourv fc Wickramasinghd 



2000[ INelemans et al.ll2001at iWood et al.ll2002ft and in 

AM CVn population synthesis s tudies flNefemans et alJ 
I2001H iFarmer fc Phinnevl 12003]) . However, the use of 
discrete, single- valued M-R relations limits what we can 
learn about these systems. For example, they do not pro- 
vide a means of inferring internal donor properties from 
observations. Nor do they allow us to model the effects 
of the expected differences between systems formed in 
a given channel. For example, in the WD-channel there 
shou ld be variations in the donor's entropy at ini tial con- 
tact (|Tutukov fc Yungelsodll99lllBildstenll2002|K In the 
He-star channel, the amount of He burned to C impacts 
the syste m's minimum P nr h and the donor's subsequent 
evolution ijErgma fc Fedoroval Il990|) . A discrete set of 
M-R relations does not allow us to model the impact 
of these differences on the AM CVn population nor to 
interpret potential observational signs of such variations. 

Progress on these fronts has been hampered by a lack 
of models in the mass range (M2 < 0.1M Q ) applicable 
to the AM CVn donors. Our recent calculation of mod- 
els for lo w-mass donors in ultraco mpact low-mass X-ray 
binaries l|Delove fc Bildsten! 120031 hereafter DB03) now 
provides a set of models continuously parameterized by 
the donor's specific entropy in this regime. In contrast to 
the above M-R relations, our model set's M-R relation is 
not single- valued, therefore allowing us to make connec- 
tions between the macroscopic and internal properties of 
the donor. In this paper we apply the DB03 models to 
the AM CVn binaries. 



In §2, we summarize the details of our models and com- 
pare their continuous M-R relations to the discrete set 
of M-R relations above. We highlight how a range in 
the donor's physical parameters translates into a range 
in observables, focusing on the relation between the sec- 
ular mass transfer rate, M, and P or b- We apply our 
models to the known AM CVn binaries in §3. For sys- 
tems where q = M2/-M1 is known or inferred, there is 
a single- valued relation between M and M^. Making 
reasonable assumptions about the nature of the donor 
and accretor, we place limits on M and Mi. Our mod- 
els provide a connection between these constraints and 
limits on the donor's specific entropy, the first such con- 
straints on the AM CVn donors. With the current ob- 
servations, this does not yet provide a strong constraint 
on donor properties (except possibly for the case of V803 
Cen which must have T c < 10 6 K given its assumed q). 
With stronger observational limits, we may be able to 
constrain the donor's prior evolution and possibly pro- 
vide a means of determining individual system's forma- 
tion channel. In §4, we use our models to explore how a 
range in donor specific entropy at contact in WD-channel 
AM CVn systems alters the M, P or b distribution from 
that calculated assuming a T = M-R relation. A siz- 
able fraction of systems in this population make con- 
tact withi n several 100 Myr or so of t he end of the last 
CE phase l|Tutukov fc Yungelsonl 1 9*961 . not allowing the 
donors to cool much before mass transfer begins. These 
high specific entropy donors follow tracks that differ from 
the T = evolution track. The consequence is a popula- 
tion that occupies a larger region of M-P or b parameter 
space. In general, higher entropy donors follow tracks 
shifted towards longer P or b relative to fully degenerate 
donors. We also discuss how this calculation impacts 
the gravity wave, GW, signal produced by this popula- 
tion. In particular, by the shifting of systems with hot- 
ter donors to longer P or b, the number of AM CVn sys- 
tems, relative to a population with only fully degenerate 
donors, that will be undetectable behind the detached 
WD- WD binary GW background increases. Finally, in 
§5 we summarize our results and discuss future directions 
for study. 

2. HOT LOW-MASS HELIUM DONORS: IMPLICATIONS 
FOR INTERACTING BINARIES 

Details concerning the ca lculation of our low-mas s 
He models can be found in iDelove fc Bildsten! IpTOl . 
Briefly, we integrate hydros tatic balance using a mod- 
ern equation of state (EOS ) ijChabrier fc PotekhInlll99St 
iPotekhin fc Chabrierll200"oT) that includes realistic treat- 
ment of Coulomb contributions, which are important at 
the low densities of these low-mass objects. Knowing the 
internal thermal profile of an object is impossible without 
following a given evolutionary scenario, so we assume the 
limiting case of an adiabatic internal profile; sec DB03 
for further justification. We calculate models for 2 < 
log(T c /K) < 7.9 and 2.0 - 2.3 < log(p c /gcm" 3 ) < 6.6 
(the lower limit depending on T c ) to produce a set of 
models parameterized by specific entropy. 

We show in Figurc^ne M-R relations for our He mod- 
els by the solid lines. Each line shows a sequence with 
the same T c and terminates with models that fill their 
Rl at P or b « 170 min. The coldest isotherm (T c = 10 4 
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Fig. 1. — A comparison between the He WD models of DB03 and 
the M-R, relations currently used in the AM CVn literature. The 
solid lines show a set of isotherms for the DB03 He WD models 
with log(T c /K) between 4.0 and 7.5 at inte rvals of Alog(T c /K) 
=0.5. The lower dashed line shows the T = O lZapolskv fc Salpeterl 
tlMflh He WD relation. The upper das hed line shows the M-R 
relations for the He-star channel used by [Nclemaiis ct alj 12001bl) . 
The dotted lines show the R^ filling solutions at the P or b of the 
indicated systems (see Table HI for these system's parameter). 



K) is equivalent to the T — M-R relation. The differ- 
ence in this curve and the iZapolskv fe Salpeterl l|1969D 
relation as fit by lNelemanset*luTl|200Tb(i (the lower 
dashed line) is due to the more simplified treatm ent of 



Coulomb physics in the 
as compared to that o f 



Zapolskv fe Salpeterl lll96l EOS 
Chabrier fc Potekhinl l)1998|) and 



iPotekhin fc Chabrierl l|2000D . The turn-over in the M-R 
relations seen at low T c is due to Coulomb interactions 
becoming the dominant source of binding in the WD. 
As T c increases, thermal contributions to the EOS at 
low density eventually dominate Coulomb contributions 
and the M-R relation no longer turns over (for He ob- 
jects, this occurs for T c « 10 5 K, DB03). For higher T c , 
there is a mini mum mass below which equilibrium mod- 
els do not exist llCox fc Guililll968t ICox fc Salpetelll96l 
lHansen fc Soangenberd 119711 iRappaport fc Jossl Il984fl 
that is related to the transition between thermal pres- 
sure and degeneracy pressure dom inating the object's 
structure l|Delove fc Bildstenl 12003]) . Correspondingly, 
for each isotherm, there are 2 branches to the M-R re- 
lation: an upper branch where the object is thermally 
supported and a lower branch where degenerate elec- 
trons provide the dominant pressure support. For a 
fixed M2, on the lower branch, R2 increases with T c , 
while on the upper branch, R2 decreases with increas- 
ing T c (a fact related to the negative heat capacity of 
thermally supported stars). For comparison, the upper 
dashed line shows the iNejgman^t^y (|2QQlb l) fit to the 
He-st ar M-R track oflTiitukov fc Fedoroval l|1989T) used 
in the lNelemans et alJ (|2001b|) population synthesis cal- 
culation and elsewhere. 

The dotted lines in Figure ^ show M-R relations re- 
quired for the donor to fill it Roche lobe (equation |2J) at 
the P orD for the indicated AM CVn systems. At fixed 
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Fig. 2. — The relation between M and P OT h for fixed Mi = 
0.6Mq as produced by our He WD models. The solid lines show 
this relation at constant T c = 10* , 10 6 , 5 X 10 6 , and 10 7 K (from 
bottom to top). The symbols indicate where along each isotherm 
M 2 = O.O1M (triangles), 0.02M Q (squares), and 0.05M Q (pen- 
tagons). Along several hotter isotherms, the same mass is indicated 
at two points, a situation resulting from the existence of both up- 
per and lower branch models at these T c . In these cases, for the 
shorter P or b point, the donor is on the lower M-R branch, while for 
the longer P nv h po i nt, th e donor is on the upper M-R branch. The 
Tsugawa & Osaki ( 1997) He disk stability criteria (with q = 0.05) 
are shown by the dashed lines. At longer P or b, allowing for hot 
donors produces a several orders of magnitude range in M. 



P or b, the minimum M2 a donor can have and satisfy this 
constraint is set by the intersection of each dotted line 
with the T = M-R relation. For M2 greater than this 
minimum, there is a continuous set of Rl filling solu- 
tions parameterized by the donor's entropy. Thus, at 
fixed P rb; the donor's specific entropy determines M2. 
On the lower branch of the M-R relations, hotter donors 
must be more massive as seen along the dotted lines in 
Figure n] 

The fact that M2 and R2 can take on a range of values 
produces a range in binary parameters that is related to 
the donor's specific entropy. To illustrate this, consider 
M, which, assuming conservative mass transfer, is given 

by' 



M 

~M~2 



1 



+ 5/3-2g 



(3) 



Here J is the orbital angular momentum, J < 
the angular momentum loss rate from GW emission 
ijLandau fc Li fshitz 1962) and 



d lnP 2 
<91nM 2 



(4) 



Throughout this paper we assume that the donor re- 
sponds adiabatically so that n is calculated along a well- 
defined adiabat. This is equivalent to assuming that no 
specific entropy is generated or lost in the donor and that 
M is high enough that the orbital evolution timescale is 
shorter than the donor's cooling time. This latter as- 
sumption seem s well justified for ultracompact binaries 
l)Bildstenll2002T) . but see §4 for a further discussion of the 
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Fig. 3. — The relation between M and P or b along evolution- 
ary tracks. Each line shows evolution from an initial configu- 
ration of M 2 ,i = O.2M and Mi ;i = O.6M . The solid lines 
show the adiabatic evolution of our He WD models initially with 
T c = 10 4 , 10 7 , 3 X 10 7 , and 6 X 10 7 K (left to right) on the lower 
branch of the M-R relations. The dot-dashed lines show the same 
for initial T c = 6 X 10 7 and 5 X 10 7 K (a gain, left to right) on 
the upper branch. The dashed lines show the evolution using the 
INelemans et alj 12001b!) fits for the T = (lower) and He-star (up- 
per) M-R relations. 



validity of this assumption. With equations J2J) and ©, 
we relate M to P or b based on M 2 , R2, n, and assuming 
some value for Mi. 

First consider this relation assuming there is no restric- 
tion on T c (other than it is low enough so that He burning 
cannot occur) at any P or b- In this case, M can vary by 
as much a factor of 100 or more at fixed P or b- This is 
illustrated in Figure [3 by the solid lines, which show the 
M-P or b relation along lines of constant donor T c assum- 
ing a fixed Mi = O.6A/0 (M varies a small amount for 
Mi = 0.4 — IAMq). It should be noted that these lines 
do not represent an evolutionary sequence, since n along 
an adiabat is used to calculate M. Instead, they show 
the parameters needed for a system to have some M at 
a given P or b- The dashed lines show the upper and lower 
thermal stability criteria for a He disk as determined by 
iTsueawa fc Osakil i|1997|) . Above the upper line, the disk 
is in a high-luminosity stable state; below the lower line, 
the disk is in a low-luminosity stable state. For systems 
located between these two lines, the disk is expected to 
exhibit periodic outbursts (see the discussion in §3). A 
population of AM CVn binaries with variable evolution 
histories should be expected to have a range in M at 
fixed Porb, with M providing a diagnostic of M2 and T c 
in individual systems. Additionally, the rate at which a 
given AM CVn systems evolves in M2 and in P or b also 
varies with the donor's specific entropy (DB03). 

The plausible range in T c for a given P or b depends 
on the formation channel and past evolution so that the 
range in M at longer P or b will depend on the properties 
of the progenitor population. We show this by consid- 
ering adiabatic evolution from the minimum P or b- The 



donor in this case follows an adiabat set by its specific 
entropy at contact. The core temperature decreases with 
M2 as the donor expands along the adiabat. We show 
a few of the resulting M-P or b evolution tracks in Figure 
E3 calculated from the initial conditions M± t i = O.6M0 
M.2,\ = Q.2Mq. The solid lines show tracks for our He 
models starting on the lower, degenerate, branch of the 
M-R relations, each line corresponding to a different ini- 
tial T c (between 10 4 and 6 x 10 7 K). The dot-dashed 
show hot models initially on the upper branch. The 
dashed lines show the evolution along the He-star (up- 
per dashed line) and T = (lower dashed line) tracks of 
INelemans et alJ (|2001b|) . Donors from the WD-channel 
are much less likely to make contact while on the upper 
branch (due to their rapid cooling while hot). Objects 
that do make contact while on the upper branch will fol- 
low tracks similar to the dot-dashed lines at much higher 
M than the those with T = donors (shown by the 
lowest solid line Figure EJ. The distribution of M at a 
specific P or b depends on the distribution of the donors' 
specific entropy at contact (assuming donors evolve adi- 
abatically). We discuss in §4 the properties of this dis- 
tribution in the c ontext of the WD-chann el population 
synthesis model of INelemans et all j|2001b|) . 

The AM CVns produced through the He-star channel 
will have, on average, donors with higher specific entropy 
than those from the WD-channel, thus populating a re- 
gion in Figu re EH above the T = tr ack. The track pro- 
duced by the lNelemans et aTl l)2001b|) fit to a He-star evo- 
lution (the upper dashed line in Figure C3 provides one 
such example. The difference in the slope of this track 
and the DB03 adiabats may be due to several causes — the 
donor cooling over the course of the lTutukov fc Fedoroval 
l)1989j) calculation is one possibility, differences in the in- 
ternal structure of the models another — and indicates 
that the approximations inherent in our models and in 
calculating their evolution are only a first step towards 
a more realistic understanding of compact binary evolu- 
tion. We will discuss some future directions for improve- 
ments to our calculations below. 

3. APPLICATION TO THE KNOWN AM CVN SYSTEMS 

We now apply our models to the AM CVn systems 
and determine the constraints that observations place 
on the binary and the donor. Each binary has a range 
in Mi over which the donor can fill its Roche Lobe, 
leading to a range of possible binary parameter. For 
most AM CVn systems, q = M2/M1 has eith e r been 
spect roscopically measured (|Nather et all Il98ll iMarshl 
Il999|) or can be inferred from P or b and the so-called 
supcrhu mp period observed in many of the AM CVn 
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well as Table [Q) . We use these q to obtain a unique 
relation between M and M2 for each AM CVn system. 
These relations, calculated using equation ©, arc shown 
in Figure 0] by the solid lines. 

The minimum value of M2 for each system is set by 
the intersection of the DB03 T = M-R relation with 
the Roche-lobe filling solution to equation at the ap- 
propriate P or b (see Figure IT). The upper limit on M2 is 
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Fig. 4. — For the systems listed in Tabic Q the solid lines show 
the relation between M2 and M based on P or b and q = M2/M1. 
The lower limit on each curve is set by our T = M-R relations, 
the upper limit by the requirement that M\ < IAMq. For some 
systems, the accretion disk's thermal stability, inferred from the 
system's photometric behavior, further constrains M. The bold 
lines with arrows indicate such upper and lower limits where appli- 
cable (for clarity we show only one arrow for AM CVn and HP Lib 
which share a lower disk-set M limit). The circles along each solid 
line indicate where our donor models have T c = 10 5 , 10 6 , and 10 7 
K (from left to right). 



set by q and the requirement that Mi < IAMq, the 
Chandrasekhar mass. The relations described by the 
solid lines only depend on the donor model through n 
in equation ©. For M 2 < O.1M , DB03 show that 
n varies by w 5% between He models with T c = 
and T c = 10 7 K. Hence, these M-M 2 relations are only 
weakly dependent on the donor model. Between these 
limits, each system has a potential spread in M 2 of a 
factor of about 3. The exception is V803 Cen, whose 
very small q sets a minimum whose accretor mass of 
Mi }ln i n « 1.3-M©, so that M 2 has an extremely narrow 
allowed range: 0.021Af Q < M 2 < 0.022M Q . However, 
since there is some uncertainty in the iden tification of 
Pnrh and the superhump period in V803 Cen ijKato et alJ 
|2004|) . and therefore in the derived q, for this system, 
these constraints should be considered tentative. For the 
other systems, roughly half (HP Lib, CR Boo, CP Eri, 
and GP Com) have M\,mux > 0-5M<n , indicating that 
the a ccretor could be a C/O or Hybrid ijlben fc Tutukovl 
IT985T) WD. However, since Mi may have grown as much 
as « 0.1 - O.2M during the AM CVn phase, the pos- 
sibility the accretor in these systems is a He WD is not 
ruled out. 

Observations place further constraints on M 2 in some 
of these systems. There is an observed correlation be- 
tween the photometric behavior of AM CVn systems and 
their P rb- The short period systems, AM CVn and HP 
Lib, ar e bright and do not exhibit large amplitude vari- 
abilit y <jWarnerll995HSolheim et al.ll998tlSkillman et all 
Il999j) . The intermediate period objects — CR Boo, KL 
Dra, V803 Cen, and CP E ri — exhibit large amplitude 
variability (« 3-4 mag nitude: lWood et al.ll987tlSchwartd 



19981 IWood et all 120021 IQ'Donoghue fc Kilkenny! fl989t 
Patterson et all 120001: lAbbott et alJ 11992(1 . while the 
long period systems, GP Com and CE 315, are con- 
sistently faint and again do not exhibit large ampli- 
tude brightening jRicher et al.l fl973t IRuiz et all 120011: 
IWoudt fc WarneJl2001(l . This behavior is believed to be 
explained with a thermal-tidal He accretion disk instabil- 
ity model l |SmaUH983t iTsugawa fc Osakilll997(l that re- 
lates M to these different obs e rved p hotometric regimes. 
Taking the lTsueawa fc Osakil l|1997(l results, we use their 
stability criteria to place limits on the allowed M in the 
AM CVn systems based on the photometric behavior of 
each. In some cases, disk stability gives stronger con- 
straints than those above and we show these additional 
limits in Figure 0] as upper and lower limits on M. 

From our He WD models, we can also determine T c as 
a function of M 2 in each system. We indicate the models 
with T c = 10 5 , 10 6 , and 10 7 K by open circles along each 
system's M-M 2 relation in Figure 0] We summarize the 
limits placed on these systems from donor and accretor 
structural constraints, along with other system proper- 
ties, in Table ^ For those systems in which He disk 
stability arguments provide stronger limits, we provide 
these limits in Table |5J These represent the first limits 
placed on the donor's specific entropy in AM CVn sys- 
tems. Again, given current observations, most of these 
limits are not strong. However, the lower limits on AM 
CVn and HP Lib set by disk stability arguments indicate 
their donors are far from being T = objects. In V803 
Cen, the donor has a maximum T c w 10 6 K, placing it 
on the degenerate branch of the M-R relations. Perhaps 
this provides evidence that V803 Cen formed through 
the WD-channel; however, a better understanding of the 
range in the donor's initial specific entropy in each for- 
mation channel is needed to support such a claim. Com- 
bined with stronger observational limits on other AM 
CVn systems, connections between the formation chan- 
nel and the donor's specific entropy could provide one 
tool for determining a system's prior evolution. 

4. AM CVN SYSTEMS FORMED THROUGH THE WD 
CHANNEL 

We have shown that the AM CVn binaries can oc- 
cupy a larger region of M-P or b parameter space than 
previously thought. There remains the question, how- 
ever, of how much of this parameter space they ac- 
tually occupy given reasonable evolution scenarios. In 
the WD channel, variations in M 2 and in the time the 
donor has to cool before it makes contact, t cont , produces 
donors with different initial specific entropies. Figure [21 
showed how this can lead to a wide range of evolution- 
ary tracks. In the He-star channel, the donor initiates 
mass transfer as a He-burning star and its structure de- 
pends on M 2 and the amount of He burned to C/O. If 
the composition was fixed, then there would be a sin- 
gle M-R relation (the He-burning main-sequence) and 
all systems would have the same minimum P or b and the 
same subsequent evolution, modulo the variations due 
to a range in Mi (up to a factor of 5). However, vari- 
ations in t cont impact how much He is processed before 
nuclear reactions are quenched by mass loss, affecting 
the donor's M-R relation, the binary's minimum P or b, 
and producing variation in the subsequent evolution dur- 
ing the AM CVn phase. Our current models cannot ac- 
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Table 1. AM CVn System Parameters with Limits Set by Donor and 
Accretor Stellar Structure 



System P orb (s) P ah (s) 



'M 2 ,, 



log | 



M 



log I 



References 



AM CVn 


1029 


1051 


0.087 


0.035 


0.40 


-9.43 


- -8.01 


7.48 


HP Lib 


1103 


1119 


0.056 


0.032 


0.625 


-9.49 


- -8.46 


7.28 


CR Boo 


1471 


1487 


0.042 


0.023 


0.547 


-10.12 


- -9.09 


7.10 


KL Dra 


1502 


1533 


0.081 


0.022 


0.277 


-10.37 


- -8.47 


7.40 


V803 Cen 


1612 


1619 


0.016 


0.021 


1.29 


-10.12 


- -10.02 


5.95 


CP Eri 


1701 


1716 


0.03 


0.019 


0.641 


-10.46 


- -9.46 


6.98 


CP Com 


2794 




0.022 b 


0.010 


0.46 


-11.66 


- -10.50 


6.70 


CE-315 


3906 ± 42 




0.0125 b 


0.006 


0.48 


-12.47 


- -11.28 


6.48 



Nclcmans^^aL^ (2001a) 
Skillman et al. (1999 , ) 
Solheim et*a!*"(1998) 
Patterson^ral"75u02 ) 
Provencal et 
^ooTTTt^u" 




Patterson et 
et al. 
Pattersoije^al. 
Patterson f 2001V " 
Marsh (1999) 
Ruiz et al. (2001) 
D. Steeghs, priv. comm. 



a Lower limits on Ma, Mi, M, and T c 
b Determined from spectroscopy. 



are set by DB03 T = He M-R relations, upper limits by Mi < 1.4M . 



Table 2. Limits on AM CVn System Parameters set by He Disk 
Stability 21 



System (^) log (_^ =r ) log(^) 

AM CVn 0.039 0.45 -9.30 - -8.01 b 6.38 - 7.48 b 

HP Lib 0.038 0.677 -9.30- -8.46 b 6.52 - 7.28 b 

KL Dra 0.023 0.282 -10.34 - -9.00 5.33 -7.22 

GP Com ■•■ •■■ -11. 66 b - -11.14 < 6.32 



a Calculated using Tsueawa & Osaki ( 1997) He disk thermal stability criteria. 
b Limit set by donor or accretor structural limits; see Table HI 



commodate varying internal composition and so we will 
only address vari ability within the WD-channel. W e re- 
fer the reader to iPodsiadlowski eta!] 1)20021 |2003|) and 
iNelson fc Rapnanortl l)2003l) for discussions of the evolu- 
tionary variations in AM CVns and other ultracompact 
binaries formed through the H-depleted main sequence 
star channel. 

We now show how a range in i con t affects the P OI b-M 
distribution of WD-channel AM CVn sy stems. We begin 
with th e population synthesis d ata from INelem ans et al] 
l|2001b() . We refer the reader to lNelemans et all (|2001bf) 
and references therein for the details of their model. The 
information needed for our calculation are the initial ac- 
cretor and donor masses, M\ t i and M2 t i, at the beginning 
of the AM CVn stage, t con t, and the galaxy's age when 
the system begins the AM CVn stage. We use i con t to 
determine the donor's specific entropy at contact using 
existing He WD evolution calculations. We evolve each 
system adiabatically from contact assuming conservative 
mass transfer driven by GW emission up to the present 
age of the galaxy (taken to be 10 Gyrs) in agreement with 
iNelemans et al.l (|2001bfl . We removed systems where Mi 
grew larger than 1.4 Mq during this evolution. 

4.1. Determining the Initial Donor Model 

The donor starts out as the core of a red giant branch 
star which initiates unstable mass transfer and a CE 
phase. The state of the outer regions of the donor af- 
ter the CE phase ends (e.g., how much of a H envelope is 
left) will depend on the uncertain details of CE evolution, 



but the central conditions of the donor will be largely un- 
changed from their state at the start of the CE phase. 
Afterwards, the donor will settle towards a cooling He 
WD configuration on its thermal time (~ 10 6 -10 7 yr). 
In the meantime, it cools and contracts until it makes 
contact at a time t caat later. The system's evolution 
through these phases determines the specific entropy of 
the donor at contact, and the distribution of the donor's 
specific entropy at contact determines the number of sys- 
tems evolving along a specified range of adiabatic tracks 
(see Figure |2J during the subsequent AM CVn phase 
and, hence, the resulting population distribution. For 
donors with M2 « 0.01 — 0.03M Q , as is appropriate for 
the known AM CVns, it is the specific entropy of the 
inner ss 10% of the initial donor mass that matters for 
current the evolution. We thus emphasize a proper ac- 
counting of the donor's central conditions throughout the 
pre-contact evolution to ensure we calculate a realistic 
population distribution for P rt> ^ 20 min. 

To determine the initial donor model, we begin by de- 
termining the central conditions in the d onor's progenitor 
at the time the CE phase begins. The INelemans et alJ 
<|2001b|) data gives the progenitor's age and He core mass 
at this time. These two quantities constrain the 
progenitor's initial mass. We use the EZ stellar evolu- 
tion code (der ived from Pet er Eggleton's stellar evolu- 
tion code; see lPaxtonll2004 and references therein) to 
calculate p c , T c in a star of this initial mass and age. 
The central conditions in the donor do not change dur- 
ing the rapid CE phase and we use the p c , T c just found 
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to determine the initial conditions for the WD cooling 
phase. To calculate the WD cooling, we use the pure- 
He WD cooling models of lAlthaus fc Benvenutol (|1997fl 
(hereafter AB97). For a given WD mass, these models 
do not necessarily have a (p c , T c ) pair matching that 
just found for the core (due to differences in the interior 
profile between the two models). We map from the core 
model to the WD cooling model by fixing p c and M 2j i, 
which then determines T c for the cooling model. While 
this is an ad hoc methodology, without being able to cal- 
culate the thermal evolution of the donor relaxing from 
its post-CE configuration to the WD configuration it is 
difficult to argue for a more accurate prescription. In 
any case, the change in central specific entropy involved 
in this mapping is usually small compared to the en- 
tropy lost during the cooling. We then use t cont and the 
AB97 cooling sequences to determine the donor's central 
conditions at contact. Finally, to determine the start- 
ing AM CVn evolution model, we find the DB03 model 
with Mi = M2 t i and a specific entropy that equals the 
central specific entropy of the AB97 model at contact. 
Generally speaking, most donors make contact as mildly 
to extremely degenerate objects, but there is a fraction 
of non-degenerate donors in the population. Once the 
initial model is determined, we calculate a from the re- 
quirement that Rl — i?2 at contact and integrate the 
system's evolution with M calculated from equation © . 
We calculate P or b from equation (0) and n from our mod- 
els assuming adiabatic evolution. 

While providing a reasonable estimate of the central 
conditions in the donor's core, this prescription does not 
allow us to determine the evolution of the donor's outer 
regions prior to contact. This will impact the evolution 
of these AM CVn systems at short P or b- For example, 
i?2 (at fixed T c and M2) is smaller in our iscntropic mod- 
els than in the thermal equilibrium AB97 models, as is 
expected. If we were to take the AB97 model's R2 at 
contact instead of the DB03 R2, the initial P or b would 
increase by up to a factor of 2-3, depending on the donor's 
degeneracy. Additionally, using iscntropic donors post- 
contact produces a different evolution than obtained by 
self-consistently evolving the structure of the donor at 
contact (something we are not able to do with our cur- 
rent models) . The discrepancy between the two evolution 
tracks should grow smaller as the donor loses mass and 
only the central regions of the star matter. 

Other factors also create uncertainty in the evolution 
during the early AM CVn phase. The cooling rate of a 
He WD depends on the amount of H left in its envelope 
after the CE event a nd whether H/He d i ffusive mixing 
is considered or not jPrieb e et al.||I999^ |Althaus et al" 



2001f|Althaus fc Benvenutolll997t iBenvenuto fc Althaui 
19981 iSarna et al.ll2000ft .' The amount of H that remains 
after the CE is uncertain (see, e.g.. llben fc Livioill993D . 
In general, more H leads to slower WD cooling due to 
its greater opacity and the possible residual H burning 
in the lower atmosphere. Further uncertainty is added 
by t he questions of whether strong shell flashes can oc- 
cur llDriebe et al.lll999t ISarna et al.ll2000t lAlthaus et all 



12001 ) , altering the subsequent cooling rates by rapidly 
consuming large amounts of H. All of this leads to signif- 
icant uncertainty in the AM CVn population distribution 
at short P or b- However, pure He WDs cool the fastest 



and the isentropic donors have the most compact con- 
figuration for a given central specific entropy. Therefore 
our results provide a lower limit on the impact of finite 
T c on the WD-channel AM CVn population. 

The uncertainties associated with the CE phase and 
the subsequent He WD cooling phase lead to uncertain- 
ties in i?2 near contact. While this clearly impacts the 
short period AM CVn distribution, it may also be signif- 
icant to the overall population. Donors with larger R2 at 
contact — as expected for hotter, isothermal donors — are 
not as likely to experience a mass transfer instability due 
to advection of angular momentum onto the accretor (ei- 
ther during a direct-impact accretion phase or because 
the accretion disk at short periods has a small radial dy - 
namic range : INelemans et al.ll2001bl iMarsh et af1l2004() . 
Taking account of this effect could, therefore, increase 
the number of WD-channel systems that survive to be- 
come AM CVn binaries (as well as the upper limit on 
M 2 in those that do). 

4.2. The Resulting Population 

We now compare the WD-channcl population calcu- 
lated using the DB03 T = M-R relations with the 
one calculated taking into account the pre-contact donor 
cooling. We will refer to these two populations as the 
T — and RWDC (realistic WD cooling) populations, 
respectively. This T = population differs slig htly from 
the analogous results of Ncl emans et alJ l(2001bD . as their 
fully degenerate M-R relation differs from DB03. In 
the M2 range of interest, the iNelemans et alJ (|2001bf) 
donors are more compact than DB03's so that the re- 
sulting M(Porb) relation from the DB03 donors is shifted 
upward by ~ 10 — 20% from that for the INelemans et all 
l)2001bD donors. In the RWDC population, the track each 
system follows is determined by the donor's specific en- 
tropy at contact, with the relative number of systems 
evolving along each a diabat determined by th e distri- 
bution in tcont- In the INelemans et all (|2001bD data w 
10%, 30%, and 70% of the systems have t cont < 100 Myr, 
500 Myr, and 1 Gyr, respectively, producing a significant 
fraction of systems not evolving along the T = track. 
Since these numbers are based on using T = radii for 
the pre-contact donor, they are a lower limits to a self- 
consistcntly determined distribution of i C ont, where more 
systems would have hotter donors at contact. The distri- 
bution of tcont also depends on the assumptions made in 
calculating the CE phase. A larger CE efficiency parame- 
ter (the ratio of the change in envelope b inding energy to 
changes in orbital energy) than used bv INelemans et all 
l)2001b|) . for example, would produce systems with wider 
separations at the end of the CE phase and reduce the 
number of systems with hot donors. 

Figure [S] shows the distribution of AM CVn binaries 
versus P orD in the two populations. Both exhibit, at first, 
the rapid increase in the number of systems with increas- 
ing P or b due to the increase in binary evolution time with 
P or b- The T = population (dot-dashed lines in Figure 
has a sharp cut-off at P or b ~ 80 min set by how far 
a T = donor can evolve in 10 Gyr. Hot donors make 
contact at, and evolve out to, longer P or b than fully de- 
generate ones. The range in the donor's specific entropy 
in the RWDC population therefore prevents a sharp cut- 
off in systems above P or b ~ 80 min. Instead there is a 
factor of s=s 20 decrease between the number of systems 
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100 



c 




4.5 



P orb (min) 

Fig. 5. — A histogram of the number of systems as a function of 
P or b between the T = (dot-dashed line) and RWDC (solid lines) 
populations. The higher specific entropy of some donors in the 
RWDC population leaks systems towards longer orbital periods as 
compared to the T = population. The smooth curves show the 
percentage of each population lying above a given P or b- In the 
RWDC population, 65% of systems have P or b > 60 min (50% of 
the T = population lies above this period); 10% of the RWDC 
population lies above the T = population's P or b ~ 80 min cutoff. 




P„u«65 min 



-9.6 -9.4 -9.2 -9 -12.6-12.4-12.2 -12 
log(M/M yr- 1 ) 



Fig. 6. — The distribution of systems in M at the approxi- 
mate P orb of AM CVn and CE-315. Our RWDC population is 
shown by the solid li ne, the T = population is shown by the dot- 
dashed line, and the Nelemans et al. (2001b) T = result s by the 
dotted lines. A comparison between the [Nelemans c t alJ J2001bl) 
population and our T = population shows the magnitude of the 
shift in the M distribution produced by the difference between the 
Zapolskv & Saloetcr ( 1969) M-R relation and our, more accurate, 
T = relations. The existence of hotter donors in the RWDC 
produces the tail seen in that population's M distribution. At 
the longer P or b of CE-315, the spread in M increases as specific 
entropy differences have more a significant impact on lower-mass 
donor's structure. 



at Porb = 80 and the number at P or b = 100 min in the 
RWDC population. The overall result is a shifting of 
systems from shorter to longer P orD in the RWDC pop- 
ulation as compared to the T = population, resulting 
in a reduction in the number of systems we expect to 
see in the P or b range of the known AM CVn population 
(i.e. P or b < 65 min). In the RWDC population, 65 % 
of systems have P or b > 60 min, while 50% of the T = 
population lies above this period; 10% of the RWDC pop- 
ulation lies above the T = population's P or b = 80 min 
cut-off. 

The range in M that AM CVn binaries can attain at 
fixed P rb is increased in the RWDC population. Fig- 
ure ^ shows that at fixed P or b , hotter donors must be 
more massive. As M oc Mj M| (at fixed P or b and 
when Mi <C Mi), hotter donors lead to a larger M (see 
also Figure |2J) . We show in Figure [5] the distribution 

of the AM CVn WD-channel population M's at the or- 
bital periods of AM CVn, and CE-315 in the T = and 
RWDC populations. In addition, we al s o show the re- 
sults obtained with the INelemans et all l|2001b|) T = 
M-R relation to indicate the magnitude of the shift in 
M produced when using our more accurate T = re- 
lation. The M range in the T = populations is due 
to the range in M\ at each P or b- In the RWDC pop- 
ulation, the existence of higher entropy donors leads to 
systems with larger Mi at a given P or b an d the tail to- 
wards higher M seen in this distribution. The offset in 
the minimum M seen between the T = and the RWDC 
population results from the fact that for the lowest accre- 
tor masses in the RWDC population, none have donor's 
that are close to fully degenerate. It is unclear if this is 
a robust result or simply a consequence of the particu- 



lar population synthesis calculation used here to provide 
the initial conditions. The number of systems in the tail 
of the distribution increases with P or b, as a higher spe- 
cific entropy has a more significant impact on the donor's 
structure at lower Mi- The majority of systems in the 
RWDC population lie near the T = tracks as is ex- 
pected, but ~ 10% of this population have donors that 
produce mass transfer rates greater than the spread at- 
tributable to Mi variations alone (at least at longer P or b; 
see the right hand panel in Figure in particular). 

4.3. Implications for the AM CVn Gravity Wave Signal 

The galactic population of AM CVn binaries con- 
tains objects that can be detected with the planned 
Laser Interferometer Space Antenn a (LIS A) gra vity- 
wave detector mission (see |htt p: / /lisa.nasa.gov/| and 
|http://sci.esa.int/home/lisa/| tor mission details), and 
the gravity wave signal of the AM CVn binary pop- 
ulati on has been considered in several recent stud- 
ies (INelemans et alJ EoOld iFarmer fc Phinnevl 1200.1 
INelemans et al.ll2004|) . Here we discuss the differences 
between the T = and RWDC populations' gravity wave 
signals. 

A binary system in a circular orbit emits gravity waves 
at a frequency / = 2/P or b (i. e. twice the orbital f re- 
quency) and luminosity, Lqw, <|Press fc Thorneiri972l) 

32 G 4 M\Ml (Mi + Mi) . 

I^GW = ~~Z r r ; (5) 

5 c° a 

where G is the gravitational constant and c the speed of 
light. The flux, F — Lqw / '47reP , received by a detector 
a distance, d, from the source is often written in terms 
of the so-called dimensionless strain amplitude, h, given 
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Fig. 7.— An AM CVn binary's (with Afi,< = 0.6M Q , 
M2,i = 0.2Mq at 1 kpc) h versus P or b, showing the impact of 
a hot donor on the system's post-contact GW flux. The dashed 
line shows the system's pre-contact inward evolution, the sym- 
bols showing when its time to contact (for a T = donor) is 
log(t/yr) = 5.0, 5.5, 6.0, 6.5, and 7.0, left to right. The solid lines 
show the post-contact evolution for a T = donor and for donors 
hot enough to make contact at P or b = 5 (10) min. The symbols 
here show time since contact: log(t/yr) = 5.0, 5.5, 6.0, 6.5, 7.0 and 
7.5, left to right. At fixed P or b, hotter donors are more massive and 
louder GW sources. The dotted line shows the detached WD-WD 
binary confusion limit for LISA iNclcmans ct al. 20 01d) . 



by l)Press fc Thorndll97i iNelemans et al.ll2001cfl 
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(6) 

where M = (M 1 M 2 ) 3/5 /(Mi + M 2 ) 1/5 is known as the 
chirp mass. In an individual AM CVn binary, changing 
the entropy of the donor changes the relation between 
M and P or b and alters the binary's GW luminosity. We 
show a summary of how hot donors change the expected 
GW signal in Figure [7J where we consider the evolution 
of an AM CVn binary with M M = 0.6M© and M 2 ,i = 
0.2Mq located 1 kpc away. The evolution of this binary 
pre-contact is independent of the donor's state, and the 
system's ft as a function of P or b evolves inward along the 
dashed line in Figure toward shorter periods. If the 
donor is fully degenerate, the binary makes contact at 
Porb ~ 3.5 min and then evolves outward in P or b along 
the left-most solid line. The orbital period at contact 
depends on the donor's entropy. The other two solid 
lines show the evolution for systems making contact at 
Porb = 5 min (requiring a partially degenerate donor 
with T c = 5.6 x 10 7 K at contact) and 10 min (requiring a 
non-degenerate donor with T c w 6.1 x 10 7 K at contact), 
respectively. Donors with these entropies occur in the 
RWDC population (which contains systems that make 
contact at orbital periods as large as « 25 min). The 
qualitative impact of a hot donor is the increase in h as 
a function of P or b as hotter donors at fixed P or b are more 



massive, increasing M.. 

The integrated GW flux from a collection of identical 
AM CVn binaries depends on the donor's specific entropy 
in two ways. First, the flux from each system at a fixed 
/ increases with the donor's specific entropy. Second, 
each system's f = df /dt depends on n, which varies with 
both the donor's mass and specific entropy. For constant 
n in a steady-state population, one can derive a simple 
scaling for the GW energy density per logarithmic fre- 
quency interval, £ G w °c fN(f)L GW (f) cx / L GW (f)/f 
= fdEQw/df. This can be expressed in terms of the 
dimensionless characteristic amplitude in a logarithmic 
interval, h r , as 



GW 



AG 



(7) 



Taking R 2 cx M 2 , and assuming M 2 <C Mi along with 
stable mass transfer, h 2 cx f a where 

2 + 12n 

(see lPhinnevlteOOHlFarmer fc Phinnevl2003l for a deriva- 
tion of this under the assumption n = —1/3). Returning 
to Figure [7| at contact n = —0.350 for the fully de- 
generate donor, while for the donor making contact at 
5 min (10 min) , n = -0.342 (-0.333). At P orb = 20 
min (/ = 1.67 x 10~ 3 Hz), the three donors have n = 
—0.233, —0.276, and —0.308 respectively. As a is an in- 
creasing function of n and is negative for n < —0.167, 
falls off less rapidly with / for cold donors (as they 
evolve less rapidly in / than donors with higher specific 
entropy). While hotter donors produce AM CVn binaries 
that are individually more GW luminous at fixed /, their 
increased rate of / evolution somewhat mitigates this in 
the population's overall GW energy density at lower /. 

We now move from these simple analytics to a di- 
rect numeric calculation of the GW signal from both 
the T = and RWDC pop ulations. We begin, follow- 
ing INelemans et al.l l|2001c|) . by randomly distributing 
the systems in each population in galactic coordinates 
with uniform probability in galactic azimuth and with a 
probability, P{Rq, z), in galactic radius, Rg, and height 
above the mid-plane, z of 



P(R 



1 



Gi Z) 



2H 2 /3 



exp 



Rg' 


sech 




H 







(9) 



where j3 = 200 pc and H = 2.5 kpc. To calculate d 
for each system we took the position of the sun to be 
Psnn = 8.5 kp c and z mn = —30 pc, again following 
INelemans et alJ l|2001cj) . We then calculated h for each 
system from its d and current orbital parameters using 
equation ©. Figure [H] presents the resulting distribution 
of AM CVn binaries in / and h for each of these popu- 
lations. As expected, the qualitative difference between 
the T = (lower panel) and RWDC (upper panel) distri- 
butions is a shift in individual systems to lower / in the 
RWDC population with an extended tail in the RWDC 
distribution to frequencies lower than / « 4.0 x 10~ 4 Hz. 

The fraction of this AM CVn population that LISA 
will detect is set by LISA's sensitivity and the galactic 
population of detached WD-WD binaries. For a signal to 
noise (S/N) ratios of 1 and 3, LISA's predicted sensitiv- 
ity (as calculated with Shane Larson's online generation 
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Fig. 8. — The distribution in the galactic population of WD- 
channel AM CVn binaries in h and / for the RWDC (upper 
panel) and T = (lower panel) populations. The grey-scale in- 
dicates the logarithm of total number of sources in each (/, h) bin. 
The solid line shows t he WD-WD confusion limit calculated by 
INelemans et al] 1200 Id) , while the dotted lines shows LISA'S pre- 
dicte d sensitivity for a S/N r atio of 1 (lower line) and 3 (upper 
line) ILarson et al.ll20Ml200^ . 
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Fig. 9. — The distribution of LISA resolved systems in the RWDC 
(upper panel) and T = (lower panel) populations. The grey-scale 
indicates the total number of sources in each (/, h) bin. The solid 
gray lin e shows the confusion limit calculated by INelemans et alj 
12001d) . while the dotted line shows LISA'S predict ed sensitivity 
for a S/N r atio of 1 (lower line) and 3 (upper line) ILarson et alj 
12000112003) . 



script, http: / /www.srl.c altech.edu~shane/ sensitivi ty / , 
based on the results of ILarson et al. 2000, 2003) is 
shown by the dotted lines in Figure 00 Popula- 
tion synthesis studies indicate that detached WD- 
WD binaries will domin ate the GW signal at fre- 
qucncics f 5 1Q- 3 Hz llFvans Then, SmaTrl 179871 
i Lipunov et"alT 119871 iHils et all Il990t INelemans et all 
l2001d iFarmer fc Phinnevll2003|) . These systems are ex- 
pected to produce a so-called "confusion limit" , a noise 
background below 2 mHz produced by the vast num- 
ber of systems that will not be individually resolvable. 
The solid lines in Figure [3] shows the confusion limited 
noise set by this background (essentially the averaged 
GW sig nal from detached WD- WD binaries) as calcu- 
lated bv INelemans et all l)2001c|) . The WD-channel AM 
CVns in the T = population that will be detectable 
above both the LISA sensitivity curve and the confu- 
sion limited noise have an average d of 8.6 kpc. The 
corresponding RWDC population contains 13% fewer 
systems at an average d of 8.7 kpc. The decrease in 
the RWDC numbers is due mainly to the shifting of sys- 
tems towards shorter /. Even though individual RWDC 
systems are on average louder than the T = systems, 
the volume over which LISA will detect these systems 
is not significantly altered. This is due to the fact most 
of the AM CVn systems above the confusion limit in ei- 
ther population are detectable to the edge of the galaxy 
(the galactic geometry sets the lower right cutoff of both 
distributions in Figure [SJ- 

For an estimated mission length of 1 yr, LISA will be 
able to resolve frequencies separated by A/ « 1/1 yr = 
3 x 10~ 8 Hz. AM CVn binaries with an h above the 
confusion limited noise and that are the only system in a 
given resolved frequency bin will be individually resolved 
by LISA (this is a bit of an oversimplification since other 
populations, such as detached WD-WD binaries, not in- 



cluded in our calculation will also contribute to the LISA 
signal; as we are interested here in how a range in donor 
specific entropy impacts the WD-channel population, we 
ignore this added complication). Between the T = and 
RWDC populations, the distribution of LISA resolvable 
sources differs slightly, as we show in Figure [5] The main 
difference between the two distributions is the slight shift 
in the maximum of the RWDC population relative to 
that of the T = and a decrease in the total number 
of resolved systems. This is seen more clearly in Fig- 
ure 1101 were we show the histogram of resolved systems 
in each population versus /. Overall, the RWDC pop- 
ulation contains ~ 6% fewer resolved systems than the 
T = and the RWDC distribution's mean frequency is 
w 10~ 4 Hz less than the T = 0's. As compared to the 
results fo r the fully degenerate distribution calculated 
using the INelemans et al.l i|2001bf) degenerate M-R rela- 
tions, the RWDC contains k> 11% fewer systems. 

4.4. Future Calculation Refinements 

These population calculations are only a first step to- 
wards a more realistic theoretical treatment of the AM 
CVn population. Here we speculate on the possible re- 
sults of a more refined calculation. In particular, how 
would improving on the approximations necessitated by 
our current donor models alter our result? The choice 
of adiabatic evolution assumes that the donor cools on 
a time scale long compared to the mass loss time scale 
(~ M2/M). While this is almost certainly the case for 
systems near contact, it may not be the case for systems 
at long P or b . If and when adiabaticity is violated depends 
on how the donor's luminosity scales with T c and M2. If 
adiabatic evolution doesn't hold, then some fraction of 
the long P or b systems in the RWDC population could 
drop out of the AM CVn population as the donor cools 
and contracts within its Roche Lobe, reducing the num- 
ber of expected long period AM CVn systems. Whether 
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Fig. 10. — A histogram of LISA resolved sources (assuming a S/N 
of 1) versus / for the T = population (dot-dashed line) and the 
RWDC population (solid line). For comparison, the distribution 
obtaine d for the T = population when using the Nelcmans ct al. 
1200 lbl) degenerate M-R relation is shown by the dotted line. All 
distributions cut-off at the 2 mHz confusion noise cutoff. The 
RWDC distribution's mean frequency is shifted downward fs 10 — 4 
Hz from that of the DB03 T = distribution. The resolved RWDC 
population has « 6% fewer systems due to the increased number 
of systems that evolve below 2 mH z. The RWDC populati on has 
11% fewer resolved systems than the Nelcmans ct al. (2001b) T = 
population. 



GW losses can subsequently bring such systems back into 
contact will depend on how much the donor contracted 
to its fully degenerate configuration and on the binary's 
parameters. The expectation would be that the orbit at 
this point would evolve on time scales of several Gyr or 
longer, so it would be doubtful that systems falling out 
of contact would reestablish contact later. The details 
of such scenarios remains to be worked out with a more 
sophisticated calculation. 

5. SUMMARY AND DISCUSSION 

We have applied the DB03 low-mass, isentropic He WD 
models, originally calculated for the donors in ultracom- 
pact low-mass X-ray binaries to the AM CVn popula- 
tion. This model set provides a continuous relationship 
between a model's M 2 , P 2 and its specific entropy, al- 
lowing us to consider the impact of donor specific en- 
tropy on the evolution of interacting binaries. Higher 
entropy donors have larger R 2 and, for a given P or b, are 
more massive than T = objects, producing a higher 
M. Across a population of AM CVn systems, variations 
in the donor's specific entropy produces a range of M-R 
relations, leading to a range in possible Af 2 and M at 
fixed P rb- For a reasonable range in T c , this range in M 
can be up to an order of magnitude or so (Figure [SJ. 

This also means that there is a range in parameters 
allowed for the known AM CVn systems. Most systems 
have a measured (or inferred) q, which, with the mea- 
sured P or b, we used to determine a unique M 2 -M rela- 
tion for each system. Structural arguments (the donor 
cannot be more compact than a T = object nor the 
accretor more massive than 1.4 M@) and the photo- 



metric behavior of each system allow us to restrict the 
range in M%. The limits we derived allow most sys- 
tems a factor ~ 3 range in M. One possible exception 
is V803 Cen, which, with its extremely low, but uncer- 
tain q ijKato et alJ I2004T) . is indicated to have a donor 
of Mi ps 0.02Mq and an accretor with a minimum mass 
of 1.3 M ; if this is truly the case, its M ps 10 _10 Af Q 
yr _1 . Our donor models also allow us to constrain each 
donor's T c range and correlate T c with M2, providing the 
first means of inferring internal donor properties in the 
AM CVn systems. A summary of the these relations be- 
tween M2, M, and T c are shown in Figure 01 and listed 
in Tables IT] and El 

Across the galactic AM CVn population, variations in 
individual system's evolution will leave an imprint on 
the overall population distribution. Our model set is 
ideally suited to explore how variations in the time be- 
tween leaving the CE phase and initiating mass transfer 
in the AM CVn phase impacts the M-P or b distribution 
of the WD-channel AM CVn systems. Starting with 
the dat a for the WD-channel from the iNelemans et alJ 
l)20f)lb() population synthesis model, we determined each 
donor's cent ral specific en tropy using the EZ stellar evo- 
lution code i|Paxtonl|2004l) to determine the donor's ini- 
tial conditions at the end of the CE phase. We then 
used each donor's t cont and the lAlthaus fc Benvenutol 
l)1997j) He WD models to determine how much each donor 
cooled by the time the AM CVn phase started. We 
evolved each system adiabatically using the DB03 mod- 
els to determine the present day orbital configuration 
for each system in the population. The results of this 
calculation are approximate for several reasons. Calcu- 
lating the CE ph ase and He WD cooling involves sev- 
eral uncertainties (iDriebe et al.lll999t iSarna et al. 200C ; 



Althaus fc Benvenutdll997HBenvenuto fc Alt1miiJll99E: 
Althaus et aLll200lD . The use of pure He models, which 
cool faster than He WD models with H present in their 
envelopes, means our results give a lower limit on the 
initial donor entropy in the WD-channel systems. 

The DB03 models have an isentropic interior profile. 
At contact, most donors should have close to an isother- 
mal core and an entropy profile that increases outward. 
At short P or b, the M-R relations of the DB03 models 
are therefore likely not a good approximation for actual 
donors; the use of isentropic models becomes increas- 
ingly better as the donor's mass is reduced and the en- 
tropy difference between its core and outer boundary is 
reduced. Therefore, the P or b distribution < 20 min cal- 
culated here is subject to uncertainty. At periods greater 
than this, M 2 < 0.1M 2ii and isentropic models provides 
a good approximation to the actual donors that becomes 
increasingly better at lower M 2 and longer P or b- 

We compare our conservative calculation that includes 
WD cooling (the RWDC population) to the same popula- 
tion evolved along the T = track. The RWDC P orb dis- 
tribution shows slightly fewer systems at short P or b com- 
pared to the T = and has an extended tail at P or b > 80 
min, the period at which the T = population cuts off 
due to the finite age of the galaxy; 65% of the RWDC 
population has P or b > 60 min, compared to 50% of the 
T = and 10% of the RWDC population has P orb > 80 
min. At fixed P or b, the M distribution of the RWDC 
population has a tail towards higher M compared to the 



12 



Dcloye, Bildsten, & Nelemans 



T = population. The fraction of systems with M higher 
than the T = population increases with P or b since a 
higher specific entropy impacts the structure of low-mass 
donors more significantly. In addition to the differences 
between t hese two populations, t here is a difference be- 
tween the lNelemans et ail l)2001bj) T = population and 
the T = calculated with the DB03 models. In the 
Mi ran ge of interest, the radii of the iNelemans et alJ 
lj2001bj) donors are less than the DB03, introducing a 

shift in the distribution of M at fixed P or b between the 
two fully degenerate populations; the DB03 T = M 
dist ributions are shi f ted a f actor of r* 1.15 — 1.5 above 
the INelemans et alJ l|2001bj) T = distributions. We 
also calculated the gravity wave (GW) signal the RWDC 
and our T = populations would produce, assuming the 
systems were distributed throu ghout the galaxy accord- 
ing to the distribution used bv INelemans et all l)2001cj) . 
The RWDC population contains 13 % fewer systems that 
LISA will be sensitive enough to detect (assuming a S /N 
of 1) with a large enough h to be seen above the unre- 
solved background of detached WD- WD binaries. This 
is due to the increased number of systems in the RWDC 
population evolving out to longer P or b and being lost un- 
der the confusion noise. In terms of systems that LISA 
will individually resolve, the RWDC population has 6% 
fewer systems than the DB03 T = po pulation and 11% 
fewer than the INelemans et al.1 l)2001btl population. 

One question posed to the theory community by the 
AM CVn population is why there are so few systems 



known. Currently, the most optimistic estimate of the 
local density of these systems is a factor of 5 below the 
most pessimistic estimates from theoreti c al populatio n 
synthesis models ijNelemans et al.ll200l5 iGrootl l2003|L 
The question is whether this is due to overall normaliza- 
tion problems in the population models, incompleteness 
of the known sample, a results of uncertain physics in- 
volved in binary evolutionary (such as the CE phase), or 
if there is unexplored physics involving the donor or the 
accretor that remove systems from the AM CVn popu- 
lation or caused them to appear differently from what 
we expect. Further theoretical investigations of this lat- 
ter possibility, along both the lines discussed above and 
others, will help answer these questions and be useful in 
constraining uncertainties in the physics governing the 
evolution of this, and other, binary stellar populations. 
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